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Abstract 

Following the Regge-Wheeler algorithm, we derive a radial equation for 
the brane-localized graviton absorbed/emitted by the (4 + n)-dimensional 
Schwarzschild black hole. Making use of this equation the absorption and 
emission spectra of the brane-localized graviton are computed numerically. 
Existence of the extra dimensions generally suppresses the absorption rate 
and enhances the emission rate as other spin cases. The appearance of the 
potential well, however, when n > \/£{£ + 1) — 2 — 1 in the effective potential 
makes the decreasing behavior of the total absorption with increasing n in 
the low-energy regime. The high-energy limit of the total absorption cross 
section seems to coincide with that of the brane-localized scalar cross section. 
The increasing rate of the graviton emission is very large compared to those 
of other brane-localized fields. This fact indicates that the graviton emission 
can be dominant one in the Hawking radiation of the higher-dimensional black 
holes when n is large. 
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The most striking result of the modern brane-world scenarios [1] is the emergence of the 
TeV-scale gravity when the extra dimensions exist. This fact opens the possibility to make 
the tiny black holes in the future high-energy colliders [2] on condition that the large or 
warped extra dimensions exist. In this reason the absorption and emission problems for the 
higher-dimensional black holes were extensively explored recently. 

The absorption and emission spectra for the brane-localized spin 0, 1/2 and 1 particles 
and the bulk scalar in the (4 -|- n)-dimensional Schwarzschild phase were numerically exam¬ 
ined in Ref [3]. The main motivation of Ref. [3] is to check the issue [4] whether or not the 
black holes radiate mainly on the brane. The result of the numerical calculation strongly 
supports the argument made by Emparan, Horowitz and Myers(EHM), i.e. the emission on 
the brane is dominant compared to that off the brane. 

In Ref. [5] different numerical method was applied to examine the absorption and emission 
spectra of the brane-localized and bulk scalars in the background of the higher-dimensional 
charged black hole. The numerical method used in Ref. [5] is an application of the quantum 
mechanical scattering method incorporated with an analytic continuation, which was used 
in Ref. [6] to compute the spectra of the massless and massive scalars in the Ad Schwarzschild 
background. The numerical result of Ref. [5] also supports the argument made by EHM if 
the number of the extra dimensions is not too large. 

However, it was argued in Ref. [7] that the argument of EHM should be examined 
carefully in the rotating black hole background because there is an another important factor 
called superradiance [8] if the black holes have the angular momenta. Especially, in Ref. 
[7] the existence of the superradiant modes was explicitly proven in the scattering of the 
bulk scalar when the spacetime background is a 5-dimensional rotating black hole with two 
angular momentum parameters. The general criteria for the existence of the superradiant 
modes were derived in Ref. [9] when the bulk scalar, electromagnetic and gravitational waves 
are absorbed by the higher-dimensional rotating black hole with arbitrarily multiple angular 
momentum parameters. 

The Hawking radiation on the brane is also examined when the bulk black hole has an 
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angular momentum [10]. Numerical calculation supports the fact that the superradiance 
modes also exist in the brane emission. Thus the brane emission and bulk emission were 
compared with each other in Ref. [11] when the scalar held interacts with a 5-dimensional 
rotating black hole. It was shown that although the superradiance modes exist in the bulk 
emission in the wide range of energy, the energy amplihcation for the bulk scalar is extremely 
small compared to the energy amplihcation for the brane scalar. As a result the ehect of 
the superradiance modes is negligible and does not change the EHM’s claim. 

There is an another factor we should take into account carefully in the Hawking radiation 
of the higher-dimensional black holes. This is an higher spin ehect like a graviton emission. 
Since the graviton is not generally localized on the brane unlike the standard model particles, 
the argument of the EHM should be carefully checked in the graviton emission. Recently, 
the bulk graviton emission was explored in Ref. [12] when the background is an higher¬ 
dimensional Schwarzschild black hole. In this paper we would like to examine the absorption 
and emission problems for the brane-localized graviton in the background of the (4 -|- n)- 
dimensional Schwarzschild black hole. 

We start with a metric induced by a (4 -|- n)-dimensional Schwarzschild spacetime [13] 

dsj^d = —h{r)dt^ + h~^{r-)dr‘^ + r‘^{d9‘^ + sin^ 6d(j)^) (1) 

where h{r) = 1 — Of course, th is a parameter which denotes the radius of the 

event horizon and the Hawking temperature is given by Th = {n + l)/47irH- To derive 
an equation which governs the gravitational fluctuation we should change the metric itself. 
Hollowing the Regge-Wheeler algorithm [14] we introduce the gravitational fluctuation by 
adding the metric component 

dP 

5s^ = [HQ{r)dtd(j) + Hi{r-)drd(j)]P'^^ sin 6—^{cos 6) (2) 

du 

to the original metric Of course, we should assume Hq, Hi « 1 for the linearizarion. 

^This metric change corresponds to the odd parity gravitational perturbation. The even parity 
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Now, we would like to discuss how to derive the fluctuation equation. Since the metric 
(1) is an usual Schwarzschild metric when n = 0, it is a vacuum solution of the Einstein 
equation. For the Ad case, therefore, the fluctuation equation can be expressed as = 0, 
where is a Ricci tensor. For n ^ 0 case, however, the metric (1) is not a vaccnm solution. 
Thus, we should assume that the metric (1) is a non-vacuum solution of the Ad Einstein 
eqnation, i.e. where and are Einstein tensor and the non-vanishing 

energy-momentnm tensor. Adding the metric pertnrbation (2) to the original metric (1) the 
Einstein eqnation wonld be changed into 5£^i, = The main problem in the derivation 

of the fluctnation eqnation is the fact that we do not know the natnre of the matter because 
the non-vanishing energy-momentnm tensor is not originated from the real matter on the 
brane bnt appears effectively in the course of the projection of the bnlk metric to the brane. 
Thus the difficulty is how to derive without knowing the real natnre of the matter. 

Firstly, we should note that the linear order pertnrbation yields only non-vanishing S£ttf„ 
5£r^ij, and 5£g^ in the gravity side. Thus it is reasonable to assnme that the non-vanishing 
components of are dT^, dT^^ and dTg^. The corresponding components of the Einstein 
eqnation are 


4 dr^ 


'\ + h 1 dh 1 d‘^h\ icn , , , 

2r^ 2r dr A dr^ I A 


^ (3) 




h 


-h 


'dHo 2 '' 

dr r 


+ 


A 


, 1 , , 1 dh 1 d^h 

- h~^{r) A - 1 - 

2r‘^ A 2r dr A dr"^ 


Hi = de 


r(j) 


iujh ^{r)H{) — h{r) 


dHi 

dr 


^Hi = AdQe 
dr 


where A = {i + 2){i — l)/2 and 

dP 

dTt^ = dQt^{r)e"'^^ sin 9^{cos 0) (4) 

Clu 

d P 

dTrfj, = d&r(i,{r)e^‘^^ sm6—-^{cos9) 

dO 

dTg^ = d&g^{r)e'‘^^ ~ 


[15] perturbation leads a more complicate radial equation and will not be discussed here. 
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The conservation of the energy-momentnm tensor yields a constraint 


\h{r) 


dr6Qr(j>~^ (h ^{r)— —h 


dr r 


5Q 


r(j) 


Aft 

“ 74 


lU 


'dh 


dF'h 


dh r d^h\ dHi 


- (^) 31 + 331? -^0 - 33 + 333? 


2r3 


dr 2 dr"^ 


Idh (dh r d‘^h\ , , , 

2* * + 2*? +''<’■) 


2r^ \dr 2 J dr 
1 dh d^h r d^h' 

2r dr dr'^ 4 dr^ 


Hi. 


( 5 ) 


It is easy to show that as expected the rhs of Eq.(5) is zero when n = 0. Another constraint 
may be derived from the fact that only two of three equations in Eq.(3) should be linearly 
independent. However, it does not yield a new constraint because it is possible to derive 
the first equation of (3) from the remaining ones if (5) is used. Final requirement is that 
the effective potential derived from Eq.(3) should coincide with the usual Regge-Wheeler 
potential when n = 0. This restriction gives a constraint 


lim 5Qts = lim 5Qrih = hm 5Qe<h = 0. (6) 

As far as we know, we cannot make any more constraint for the energy-momentum tensor 
from the general principle. Of course the constraints (5) and (6) are not sufficient to £x <50*0, 
SQr4> and completely. Thus we cannot determine the effective potential due to the lack 
of knowledge on the matter content. Thus what we can do is to choose 50^,^ consistently 
with Eq.(5) and (6). 

We will choose 6Qr<f> = <^000 = 0 with a hope that the detailed shape of the effective 
potential does not change the graviton emissivity drastically^. Eliminating Hq in Eq.(3), 
one can derive a differential equation for solely Hi in the form 


^Of course, there are many other choices different from 5Qr<f> = <5000 = 0 which satisfy Eq.(5) 
and Eq.(6) simultaneously. Eurthermore, the different choice may yield different and probably 
more complicate effective potential. The reason we choose <50^0 = <5000 = 0 is not because it is 
more physically reasonable than other choices but because it yields an effective potential via simple 
procedure. Thus our choice is valid physically when this assumption is true. 
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In Eq.(8) we introduced a parameter e. When e = 1, Eq.(8) corresponds to the gravitational 
fluctuation and e = 0 corresponds to the scalar fluctuation. Thus one can treat the scalar 
and gravitational fluctuations in an unifled way. We used e to confirm our numerical calcu¬ 
lation which will be carried out later. Putting e = 0, we checked our numerical calculation 
reproduces the absorption and emission spectra of the brane-localized scalar calculated in 
Ref. [3,5]. 

Defining again i? = T/r, Eq.(8) is transformed into the Schrodinger-like equation 

^ - ^effW] ^ = 0 (9) 

where the effective potential is 


1/ff(r) = h{r) 


-1-1) Idh 
- 1 - 


(1 - h{r)) 


It is worthwhile noting that until this stage we have not used the explicit expression 
of h{r). Thus we can use the equations derived above for the gravitational fluctuation of 
any dd-induced spherically symmetric metric. Using an explicit expression of h{r), one can 
easily show 


Vgff(r) = h{r) 


£(£ + 1 ) 


cTutT 


where an = e[(n -|- l)(n -|- 2) -|- 2] — (n -|- 1). When n = 0, the effective potential (11) reduces 
to the usual Regge-Wheeler potential. Defining a “tortoise” coordinate z dz/dr = /i“^(r), 
the wave equation (9) simply reduces to 
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FIG. 1. The z-dependence of the effective potential for the brane-localized graviton. Fig. 1(a) 
and (b) correspond to i = 2(the lowest mode) and £ = 3 respectively. For fixed £ each potential 
with n < £ has a barrier whose height increases with increasing n. When n > £, each potential 
has a barrier in the side of the asymptotic region and a well in the opposite side. The emergence 
of this potential well in the effective potential seems to affect the low-energy behavior of the cross 
sections. 


dz^ 


+ 


UJ 


hff(' 


'S = 0. 


The effective potential for the gravitational fluctuation, i. e. 
of ^ in Fig. 1 when £ = 2 (Fig. 1(a)) and £ = 3 (Fig. 1(b)) 


( 12 ) 

e = l,is plotted as a function 
For the lowest mode £ = 2 


7 












each potential with n < i has a barrier whose height increases with increasing n. Thus the 
existence of the extra dimensions should suppress the absorption rate. Each potential with^ 
n> i has a barrier in the side of the asymptotic region and a well in the opposite side. When 
n increases, the barrier and well become higher and deeper respectively. While, therefore, 
the barrier tends to decrease the absorption rate with increasing n, the well tends to enhance 
it. Thus, the absorption cross section should be determined by the competition between the 
barrier and the well. The effective potential with £ = 3 exhibits a similar behavior as Fig. 
2(b) shows. 

The absorption cross section for i = 2 and 3 is plotted in Fig. 2 when n = 0, 1, and 2. 
Fig. 2 indicates that overall, the existence of the extra dimensions suppresses the absorption 
rate as other fields. This means that the potential barrier is more important factor than the 
potential well in determining the absorption rate. However, Fig. 2(a) shows that the cross 
section for n = 2 is larger than that for n = 0 or 1 in the low-u; region. This means that the 
potential well plays more important role in this region. 

Now we would like to explain how to compute the physical quantities related to the 
absorption and emission of the brane-localized graviton. Since the numerical method used 
in this paper was extensively used in Ref. [5,6,11], we will just comment it schematically. 
Using the dimensionless parameters x = oor and xh = ojth, the radial equation (8) can be 
written as 


x^{x^+^ - + x{x^+^ 


dx"^ 


X 


rt+l'i 

H 


) 2x^^^ + (n — l)x 


,n+l 

H 


dR 

dx 


+ 


X 


2n+4 


' + l)x 


n+l f ^n+l 


+ e {(n + l)(n + 2) + 2} - a;”+^) 


(13) 


R = 0. 


Due to the real nature of Eq.(13), R* should be a solution if i? is a solution of Eq.(13). The 
Wronskian between them is 


^The exact condition for the appearance of the potential well can be derived from the effective 


potential as n > \Jlit + 1) — 2 — 1. 





FIG. 2. The w-dependence of the partial absorption cross sections for i = 2(Fig. 2(a)) and 
(. = 3(Fig. 2(b)). Overall, the existence of the extra dimensions suppresses the absorption rate as 
in the cases of other fields. However, the reversion of the magnitude in the low-energy regime of 
Fig. 2(a) seems to be due to the emergence of the potential well in the effective potential. 


n—l 


^ dx dx - Xh^^ 


(14) 


where C is an integration constant. Since x = xh is a regnlar singnlar point of the differential 
eqnation (13), we can solve Eq.(13) to obtain a solntion which is convergent in the near¬ 
horizon regime as following: 


Qn,ti.x,XH) = ^ dc_n(x -Ih)" 


( 16 ) 


Af=0 
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where A„ = —ixu/in + 1). Of course, the recursion relation can be directly derived by 
inserting Eq.(15) into the radial equation (13). The Wronskian between Q* and Q is same 
with Eq.(14) with C = In Ref. [6] the relation between and the partial 

transmission coefficient Tn/ 

U |2 (^+ 1 / 2 )^^ 

\defl\ — -2- 

Xh 

is derived by making use of the quantum mechanical scattering theory. 

The ingoing and outgoing solutions which are convergent in the asymptotic regime can 
be derived directly also as series expressions in the form 

OO 

N=0 

The Wronskian between and JF(_) is same with Eq.(14) with C = 2i. 

Since the asymptotic solution is a linear combination of the ingoing and outgoing solu¬ 
tions, one can assnme that the scattering solution which is convergent at the asymptotic 
region is 


Tln,e{x,XH) = fY\xH)d='n,e{+){x,XH) + f!;,y {xH)d^n,e{-)i^^ ^h) 


(+)/ 


(18) 


where the coefficients are called the jost fnnctions. It should be noted that if we know 
77, the jost fnnctions can be computed by the Wronskian as following: 

n+l _ ™n+l 

= ± (19) 

In Ref. [6] the relation between the jost fnnctions and the coefficient is also derived in 
the form 


din — 


£+ 1/2 

f(-) 

Jn,l 


( 20 ) 


Combining Eq.(16) and (20) the partial absorption cross section can be written as 

7r(2£ + l)r| 


a 


BR _ 


TT 


= —(2£+l)T,= 


r2 

H 


UJ 


l/+(*ff)l= 


( 21 ) 
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where the superscript denotes the brane-localized graviton. 

The emission spectrum, i.e. the energy emitted per unit time, is given by the Hawking 
formula 


YBR 


3^BR 


00 a. 


abs 


— 1 ) 


doo 


( 22 ) 


where is a total absorption cross section. Thus the absorption and emission spectra 
can be completely computed if we know the jost functions. In Ref. [5,6,11] the function TZ 
is computed numerically from the near-horizon solution (15) by increasing the convergence 
region via the analytic continuation. Then Eq.(19) enables us to compute the jost functions 
numerically. This is an overall algorithm of our numerical technique. 



Fig. 3 


FIG. 3. The w-dependence of the total absorption cross sections for the brane-localized gravi¬ 
ton. Overall, the existence of the extra dimensions suppresses the total absorption cross sections. 
The low-energy behavior of this figure also implies the crucial effect of the potential well emerged 
when n'>l. 


The total absorption cross section is plotted in Fig. 3 when n = 0, 1 and 2. As in the 
case of the scalar held absorption the existence of the extra dimensions in general suppresses 
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the absorption rate. This indicates that the potential barrier is more important factor in 
determining the absorption cross section. However, the low-cu region of Fig. 3 seems to 
indicate that the emergence of the potential well when n > i crucially changes the order of 
magnitude. The wiggly pattern indicates that each partial absorption cross section has a 
peak point at different u. The high-energy limits of the total cross sections seem to coincide 
with the scalar case, i.e. 


lim a 

UJ^OO 


BR 

abs 


( 2 ) n + l'""’* 


(23) 



Fig. 4 


FIG. 4. The w-dependence of the emission rate when n = 0, 1, and 2. As in the cases of 
other fields the existence of the extra dimensions generally enhances the emissivities. The emission 
spectrum for n = 0 cannot be seen clearly because it is too small compared to other cases. The 
drastic change of the emission rate indicates that the graviton emission can be a dominant one in 
the Hawking radiation of the higher-dimensional black holes for large n. 


The emission spectrum of the graviton is plotted in Fig. 4. As in the case of the scalar 
emission the existence of the extra dimensions enhances the emission rate. This indicates 
that the Planck factor is most important in determining the emission rate. For n = 0 the 
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emission spectrum cannot be seen clearly in Fig. 4 because it is too small compared to the 
other cases. The drastic increase of the emission rate implies that the Hawking radiation of 
the higher spin becomes more and more important for large n. 



n = 0 

n = 1 

n = 2 

graviton 

7.67 X 10-® 

8.06 X 10-^ 

6.70 X 10-3 

scalar 

2.98 X 10-^ 

2.66 X 10-3 

1.07 X 10-2 

graviton/scalar 

0.026 

0.303 

0.624 


Table I 


The total emission rate / T^^ for graviton and scalar is summarized in Table I. Generally 
the total emission rate for the graviton is smaller than that for the scalar. But the ratio 
between them tends to increase dramatically with increasing n. Thus we can conjecture 
that for large n the graviton radiation becomes dominant in the Hawking radiation of the 
higher-dimensional black holes. 

In this paper we studied the absorption and emission of the brane-localized graviton in 
the context of the brane-world scenarios. The existence of the extra dimensions in general 
tends to suppress the total absorption cross section as the cases of other helds. In the low- 
energy region, however, the order of the magnitude is reversed (see Fig. 3). This seems to 
be due to the emergence of the well in the effective potential. The high-energy limits of 
the total cross sections for the brane-localized gravitons seem to coincide with those for the 
brane-localized massless scalars. 

As in the cases of other helds the existence of the extra dimensions generally enhances 
the emission rate. This fact indicates that the Planck factor is most dominant term in 
determining the emissivities. However, the increasing rate of the brane-localized graviton 
is very large compared to those of other helds. To show this explicitly, the increasing 
rates of the brane-localized graviton and scalar is given in Table I. This rapid increase of the 
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emissivity indicates that the graviton emission should be dominant radiation in the Hawking 
radiation of the higher-dimensional black holes when n is large. 

It is straightforward to extend our calculation to the higher-dimensional charged black 
hole background by choosing 


h{r-) = 




(24) 


where r± correspond to the outer and inner horizons. Since we have not used the explicit 
expression of h{r) when deriving Eq.(8), we can use the radial equation (8) with only choosing 
different h{r). The introduction of the inner horizon parameter r_ changes the Hawking 
temperature into 


Th = 


n+1 



_ ^n+l 


(25) 


When, therefore, r_ increases, the Hawking temperature becomes lower, which should en¬ 
hance the absorptivity and reduce the emission rate. Since the introduction of the extra 
dimensions generally suppresses the absorptivity and enhances the emission rate, the exact 
absorption and emission spectra for the charged black hole case may be determined by the 
competition between r_ and n. 

Recently, in Ref. [12] the absorption rate for the bulk graviton is studied in the 
Schwarzschild background. The fluctuations for the bulk graviton have three different modes, 
i. e. scalar, vector, and tensor modes. It is of greatly interest to apply our method to compare 
the emission rate for the brane-localized graviton with that for the bulk graviton. 
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